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Abstract: We study string interactions in the fermionic formulation of the c = 1 
matrix model. We give a precise nonperturbative description of the rolling tachyon 
state in the matrix model, and discuss S-matrix elements of the c = 1 string. As a 
first step to study string interactions, we compute the interaction of two decaying 
DO-branes in terms of free fermions. This computation is compared with the string 
theory cylinder diagram using the rolling tachyon ZZ boundary states. 
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1. Introduction 

It has been known for more than ten years that the non-critical c = 1 string can 
be described by the double scaling limit of a quantum mechanical matrix model, for 
reviews see e.g. [||, H, H. In this duality the matrix model describes discretized world- 
sheets, which become continuum world-sheets in the double scaling limit. Recently, 
a new perspective on this duality has emerged which is in spirit closer the AdS / CFT 
duality and open/closed string duality. The crucial new insight obtained in Q 
and in is that the matrix degrees of freedom of the quantum mechanical matrix 
model are to be identified with the open string tachyon degrees of freedom living on 
a system of unstable DO-branes in the c = 1 string theory. The full c = 1 string 
theory is recovered by taking a double scaling limit of the matrix model. For more 
discussion see e.g. the recent papers [0, [^, 0, |T1|, [T2[. 

The c = 1 matrix quantum mechanics can be reformulated in terms of free 
fermions moving in an external potential. After taking the double scaling limit the 
potential becomes an upside-down quadratic potential, the system of free fermions 
becomes a non-relativistic 1 + 1 dimensional field theory, and the ground state is a 
Fermi sea filled with an infinite number of fermions. The distance /i of the Fermi sea 
to the top of the potential can be identified with the inverse string coupling — 0. 
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It is an interesting question to ask how the string interactions and string diagrams 
are encoded in the free fermionic theory. In principle, one expects that different 
orders in string perturbation theory can be extracted from the matrix model by 
studying its large fi expansion. The tree-level scattering of closed string tachyons 
in the two dimensional string theory were computed in the matrix model (where 
one can compute them directly or by employing the free fermionic formulation) in 
[Q, IE . These scattering amplitudes were found to be in agreement with closed string 
computations when on the matrix model side external "legpole factors" are attached. 
These factors have poles at integer values of the external momenta, corresponding 
to discrete states in the string theory. After analytic continuation of the string 
amplitudes to Minkowski signature, these factors turn into a phase, which do not 
affect the physical scattering amphtudes. 

However, the discrete states and the legpole factors play an important role in un- 
derstanding the full structure of the c = 1 string theory. As shown in gravitational 
effects are encoded in the legpole factors. 

In the c = 1 matrix model the full perturbative S-matrix was also computed 
[0, |T^. Nonperturbatively this matrix model is unstable, since only one side of 
the potential is filled, and there are no closed-string states corresponding to the 
asymptotic region on the other side of the potential. A non-perturbatively unitary 
extension of the c = 1 matrix model is obtained by filling up both sides of the 
potential, in which case it does not describe the non-critical c = 1 string but the two 
dimensional OA/OB string theory [|T^. In this unitary theory an exact S-matrix can 
be computed [1^, [1^ . 

The partition function of the two dimensional string theory and its Kontsevich 
representation was worked out in |19| . 

Although many results have already been computed on the c = 1 matrix model 
side, in the present context we have a better understanding of the direct relation 
of the matrix model to the DO-branes of the two dimensional string theory. In 
particular, processes involving D-branes on the string side can now be studied with 
matrix model methods. 

In 0, ^, ^ the decay of a single D-brane was analyzed in detail. In this picture, 
the decay of the DO-brane can be described by a single eigenvalue rolling down 
the potential. In the fermionic language one can describe this process by a single 
fermion moving from the top of the potential to asymptotic infinity, where it becomes 
relativistic. Bosonizing the relativistic fermion leads to the identification of the D- 
brane as a coherent state of tachyons in the asymptotic region. The closed string 
interactions are encoded in the interactions of the bosonized tachyon field. A general 
open-closed duality based on this description is suggested in |^ . 



In the string theory side, the decay of the unstable D-brane can be studied using a 
description involving boundary conformal theory. The relevant boundary states have 
two non-trivial pieces, one is the Liouville part which has been exactly determined 
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in [|T^, 14|, the other is the time-dependent piece which involves a timehke scalar and 
which was first studied in |jl5|. A precise determination of the time-dependent part 
of the boundary state is actually quite subtle and we will discuss it in more detail 
later in this paper. 

As a first step in understanding in more detail the relation between the string 
worldsheet description and the dual matrix model, one would like to see how the 
string interactions are reproduced in the matrix model. In particular, since the 
matrix model has a description in terms of a simple system of free fermions in a 
quadratic potential, one would like to understand if and how string interactions 
between D-branes can be computed from the free fermionic theory. 

So far, in various papers 0, ^ ||, |21| the disk one point function describing the 
decay of a single D-brane was computed and found to agree with the bosonized 
fermionic description. A further step is made in |P, where higher Qs corrections to 
the outgoing tachyon state were discussed using the S-matrix of the of the matrix 
model. 

In the present paper we further analyze the question of string interactions in 
the presence of D-branes from the free fermionic point of view. The main point is 
to understand and compute string diagrams (the disk and the cylinder), including 
possible higher order and nonperturbative corrections, in terms of the free fermionic 
formulation of the matrix model. 

In section |] and ^ we develop the systematic formalism to compute string di- 
agrams from the free fermionic point of view. The disk diagram, which has been 
computed before serves mainly as an example to develop our formalism. As de- 
scribed in section ^ the rolling tachyon can be viewed as a combination of in and out 
states. Here we use the S-matrix formalism to compute the disk amplitude and the 
string interactions between DO-branes. We conclude that the matrix model results 
are in agreement with the boundary state computations, if in the latter we take into 
account the exchange of closed string tachyons only. In particular, we do not see 
discrete modes contributing in the matrix model computation. 

To compare the matrix model formulation of the interaction with the boundary 
state formulation, in section ^ we compute the tree level interaction between two 
decaying DO-branes. The Liouville part of the cylinder diagram is already computed 
in p3] . Tensoring with a Dirichlet state in the time direction describes for example 



the cylinder diagram between two D-instantons |jTO|, |TT|. Here we consider the tree- 
level interaction between two decaying DO-branes. We will take as our boundary 
state the one that arises by an analytic continuation of the zero mode sector of the 
Euclidean boundary state, and briefly discuss possible alternatives to this. 

The comparison of the boundary state computation and matrix model results is 
further discussed in section ^. 
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2. The interaction of D-branes : A bosonic picture 



In this section, we briefly review some background material pertaining to the c = 1 
matrix model and the asymptotic bosonization picture of free fermions. Then we 
explain how to reproduce the brane "interaction" from this picture. The Hamiltonian 
of free fermions in the inverted harmonic oscillators with Fermi energy /i is given by 



H 



[2.1) 



We can introduce new left moving and right moving chiral fermions 



where 



v{x) = \/ x^ — 2[i. 



+if,r+ij^dyv{y)^^^^^^^^ (2.2) 



:2.3) 



The factors appearing in front of and \E'/j are the WKB wavefunctions of the 
inverted harmonic oscillator. From now on, we will focus on the asymptotic region 
of large negative x, or the q = — ln{—x) — > — oo region. The reason for this will 
soon become clear. At this asymptotic infinity, the Hamiltonian is reduced to the 
Hamiltonian of relativistic fermions and we have 



2 



[2A) 



Here we suppressed 1/n corrections. This can be justified for the tree level com- 
putation in string theory. Possible systematic correction to this truncation will be 
discussed in section ^. The chiral fermions can be bosonized as 



l=:exp[-;=/'(n-W 



exp 



/ (H + d,T)dq' 



[2.5) 



This scalar boson is, roughly speaking, a closed string tachyon modulo possible dis- 
crete modes. The Hamiltonian of this scalar boson can be found easily by inserting 
( p.5|) into (|2.2| ) and evaluating ( ^TTI) . It has interactions of order O(T^) which are 
again suppressed by inverse powers of fi, 

H=\j dq{Tl' + {d,T)' + 0{T')). (2.6) 

The field T{q,t) satisfies the canonical equal time commutation relations 

[T(g,t),H(g',t)]=z5(g-g')- (2-7) 
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Near asymptotic infinity, the bosonic field T can be expressed in terms of simple 
plane waves, 

-j== Tke-'\^\'+'^'i + Tle'\^\'-'^^ . (2.8) 

The commutation relations of the creation and annihilation operators of this bosonic 
field are the usual ones, 

[n,Tl,] = \k\5{k-k'). (2.9) 
Including the normal ordering, (|2.5|) can now be written as 



/27r 

WL\(lit) = e •'-°°TfcT e m . (2.10) 

V 27r 



We distinguish and 7\., where the former acts naturally on the out vacuum and 
the latter naturally on the in vacuum. The relation between and can be found 
by comparing and equating the values of conserved quantities at very early and very 
late times. The relevant conserved quantities are related by an analytic continuation 
to the generators of the ground ring (which also correspond to conserved quantities). 
The result, as taken from 0, reads 

{Tl - T,J ■ ■ ■ {Tl - nJ6i±\k^\ ■■■±\k^\- k). (2.11) 

The leading one gives the tree level interaction, and since this expression includes all 
1/n corrections, we can include higher order corrections f^. 

States that describe the excitation of a single fermion out of the Fermi sea can 
be expressed as ipR{Qyt)\f^) or '?/'|(g, t)|/i), where denotes the Fermi sea. If we 
insert the classical trajectories of x = —A cosh t or g = — log | — log cosh t in the 
single fermion state, we get, for instance, 

^t(g,t)|/i) = -=e l/i) (2.12) 



at asymptotic future infinity t — +00. This single fermion state can be interpreted 
as outgoing coherent tachyon states emitted from the decaying D-brane 0. 

The question is how we can compute the "interaction" from the free fermion the- 
ory. The strategy is as follows. Although the theory is free, there are non-vanishing 
multipoint functions of fermions, and we will interpret the two point function of 
fermions {fi\ilj'^{q,t)tp{q,t)\fi) as an "interaction amplitude" of D-branes. We should 
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note that since each fermion in the two point function is interpreted as a decay- 
ing D-brane, to get a sensible S-matrix element of two interacting D-branes the 
(g, t) coordinates should be sent to asymptotic infinity as one usually does when 
constructing an S-matrix. Of course, this interpretation looks a bit different from 
standard quantum field theory at this moment, but in the next section we will see 
that once it is assembled together that the whole picture matches perfectly well with 
standard quantum field theory. At this asymptotic region, we can think about the 
"interaction" of two Fermi particles as 



x——e J-oojk] k j-oojfcT u) 



1 r + aa dk -ik((t + q)-{t'+q')) 1 , r oo dfe p-ifc{log A-log A') 
P~ Jo k = p~ Jo k 



-e ■>0 k = e 

27r 2tt 



(2.13) 



where in the last step we inserted the classical trajectory q = — log| — log cosh t 
in the correlation function of two Fermi fields at asymptotic infinity. We have also 
neglected the overlap of the wavefunctions of the fermions with the Fermi sea in the 
computation. The form of ( 2.13|) suggests it is the ladder summation of cylinder 



diagrams. The single cylinder can be rewritten as 

Cylinder ~ j rffcrfcue-^'^ A log A' _ (2.14) 



We see that ( |2.14| ) is composed of three parts: a massless scalar particle propagator 



and two one point functions. The scalar particle is the closed string tachyon as we 
know from ( |2.5| ). We will discuss further the match in the section ^, but already at 
this level, it is apparent that even though (|2.13|) explains the closed string tachyon 
exchange perfectly, there is no signal of discrete modes. 



3. The interaction of D-branes : A fermionic picture 

We will now try to give a more precise definition of the rolling tachyon state in the 
c = 1 matrix model. Since the rolling tachyon has a very specific behavior at early 
and late times as we have considered in the previous section, we propose that in the 
Minkowski picture it is related to a combination of an in state and an out state. 
Bulk correlators in the c = 1 matrix model are evaluated between this in and out 
state. Both the in and out state have a free parameter, which defines the behavior 
at very early times ~ aexp(— X°) and the behavior at very late times ~ 6exp(X°). 
We will therefore not view the parameters a, b as defining a single state in the matrix 
model, but as defining a pair of states consisting of an in state and an out state. 
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The in state represents an incoming coherent wave of tachyons, the out state and 
outgoing coherent state of tachyons, which agrees with the intuitive picture of the 
rolhng tachyon boundary state. 

Of course, what we say here is basically the statement that a classical fermion 
trajectory is nothing but a D-brane, but so far in the literature the rolling tachyon 
has mainly been considered as single state in the matrix model, rather than as a 
combination of an in and an out state. However, if we want to compute S-matrix 
elements of the c = 1 string, we need to specify both the initial and final conditions 
of the D-brane with the rolling tachyon, and we will see below that this is indeed the 
right picture. 

The map between closed string tachyon correlation functions and c = 1 quantities 
has been discussed in e.g. [0. It involves the asymptotic bosonization described in 
the previous section, in particular the modes in ( p.9|) are the modes of the closed 
string tachyon. The same chiral fermions that appear asymptotically will also be 
used to construct the in and out states for the rolling tachyons, but whereas closed 
string tachyons correspond to fermion bilinears, for the rolling tachyon we will use 
a single fermion operator, as discussed in The relation between the exact c = 1 
fermion and the asymptotic chiral fermions involves WKB factors as in ( p. 21) that 
need to be removed in order to obtain a finite result. This removal of WKB factors 



is also exactly what is needed to obtain finite tachyon correlators ||17|| . 

The free fermion theory in the inverted harmonic oscillator has standard even 
and odd parity real delta-function normalized wavefunctions ijj^{u!,x), as described 
in [p!7| , . We want to find the wavefunctions that have the appropriate plane wave 
behavior as x ^ — oo. The relevant wavefunctions are denoted x^(c<j,x); we could 
also have considered the other regime x +oo with corresponding wave functions 
Xr{uj,x) but will not do that here. This will be relevant once we are interested in 
the sinh boundary state (in the matrix model in which both sides of the potential are 
filled), for example. The xf, ^^re linear combinations of the standard wavefunctions, 
to be precise. 



X~L x)) \ s+(a;) s- {uj) / V ^" ^) 



(3.1) 



with 



where 



s^{uj) = ^/k([J) ± -^== (3.2) 



k{uj) = Vl + e2-^ - e"", k{uj)-^ = Vl + e2-'^ + e"^. (3.3) 
From ( ^.11 ) we readily obtain 



ip-{uj,x)J 4 V-s+(^) s-{uj) J \xliuj,x) 



(3.4) 
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so that in terms of creation an annihilation operators (with standard commutation 
relations) we have 



The asymptotics of Xl as a; — >• — oo are 
with 



^(^Us xt - s^Xl) + ^(^o^i-^^xt + s Xl) 



(3.5) 



T2i 



1 



1/2- 



F{x, (j) ^ ^\x\^ -u log \x\ + ^(u) 



with 



= I + i 



r(i-H 



r(i + iu;). 

We are interested, for the in-state, in the operator 

lim ip\t ^ -\og\x\ +\og{Xin),x) 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



acting on the ground state Similarly, for the out state, we are interested in the 
operator 

lim tjj{t = log |a:| - log(Aout), a^) (3.10) 

acting on the vacuum state To obtain a finite resulting expression, we need 
to restrict to xt (^^® appropriate left/right moving part), and to remove the WKB 
piece 

(3.11) 



1 



exp(-|a;| ). 

X\ 4 



The resulting operators are most easily written in terms of the operators 



in,a) 



out, a; 



gj<I>(uj)-jwlog^ ^ ^ ^ ^ 

2(1 + ^2Tnvy/4(^<^(^<^) ~^U}{(^U})) 

^i^{uj)-iuj log^ 

2(l + e2'^'^)V4 



(3.12) 



which will play a crucial role in what follows. They describe the modes of the asymp- 
totic chiral fermions, but in terms of the exact annihilation and creation operators of 
the c = 1 theory. The factors of log ^/^t are put in by hand here in order to facilitate 
the comparison with world-sheet results later. The anti-commutation of the fermions 



are 
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^TTLJ ^2i^{uj)—iu} log fl 



e 



„1TLU p — 2i^(w)+iL0 log fl 

The factor that appears in the the commutation relations between the in and out 
modes is the reflection coefficient which is the S-matrix of the free fermion theory. 
Most of the non-trivial physics of the c — 1 theory is contained in this reflection 
coefficient. It can be rewritten as 



Ru, 



^■KUJ ^i<^{u)—iu log fj. 



^g7ra;-ia;logM Y{\ + iuj) 



g7ra;/2-ia;logMr(l/2 + ia;) (3.14) 



. 27r 

in view of the relation 

r(l/2 + ia;)r(l/2-ia;) = — ^ . (3.15) 

coshTTo; 

The factor i?^^ decays exponentially at large negative a;, which guarantees that the 
results we obtain will be finite. This exponential decay indicates that only non- 
perturbative effects will access the region where the energy is negative. 

The refiection coefficient has nonperturbative pieces and a perturbative expan- 
sion in 1///. The perturbative expansion is of the form 



(3.16) 



with an irrelevant prefactor — ie"*'*. 

Based on our previous discussion, we now propose that the operator that creates 
a D-brane with a rolling tachyon that as a; — > — oo behaves as 

T ~ Aine-* (3.17) 



IS 

I in) 



and similarly that the out operator that annihilates a D-brane with a rolling tachyon 
that as X — > -|-oo behaves as 

T ~ Aoute* (3.19) 
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IS 



(out| ^ {l^\D^^,;,^^^ = (/.I J da;e— (3.20) 

Similarly, we can make in and out states that contain multiple branes by simply 
acting with the corresponding brane creation and annihilation operators. 

Before continuing the discussion, we first test this proposal by computing the 
energy associated to the rolling tachyon. In string theory, this should be the one 
point function of the relevant ground ring operator in the rolling tachyon boundary 
state. For the present point of view the relevant calculation is easy to write down, 
since the energy of (a^) Ms a;. The exact expression for the energy therefore is given 
by computing the expectation value of cu between the in and the out state, which 
after normalization leads to the expression 



-ia;(logAi„+logAout),^, e;^^ /£(i±M 

^ . (3.21) 

fl" p-Mlog Ain+log Aout) , / 

In this expression we absorbed a factor of ^/JI in both Ain and Aout- 

Let us try to work out this expression. In the region where uj is close to /i, 
the factor e^^/VT+~e^ can be dropped as it is equal to 1 plus nonperturbative 
corrections. We are left with 

-iu;(logAi„+logAout),,. /£(|±:^ 
f^' e-iu;(log Ai„+log Aout) , / 

If we assume the rolling eigenvalue has small overlap with the Fermi sea, we can 
extend the integral to the range from — oo to oo, and since it is rapidly fluctuating 
we can use the method of stationary phase. The fluctuating determinant is sub- 
leading in the genus expansion, so we are interested in the leading stationary phase 
approximation. Using Stirling's formula, the gamma functions behave as 



~ exp(z(a;loga; - a;) + 0{l/u;)). (3.23) 

The stationary phase is at a; = AinAout which is indeed where we expect the energy 
to be located, and the expression for the energy becomes 

^ = — 77 log(exp(-iAinAout)) = AjnAout- (3.24) 

d{-l log(AinAout)) 

This is indeed the correct energy for an orbit of the type 

X{x) ~ Aine-* + Aoute*. (3.25) 
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We can easily take the fluctuation determinant into account in this discussion. This 
merely adds an extra term z/2 to the energy, from which we see that this is indeed 
a subleading effect, since Ajn and Aout both contain a factor of ^JJi that we absorbed 
at the beginning of this calculation. 

3.1 Adding tachyons 

The construction of the tachyon operators uses crucially an asymptotic bosonization 
procedure. To add tachyons to the game, we need the normal ordered products 

p{x,t) =:^\x,t)^{x,t) : (3.26) 

in the same in and out limits as above. This operator in boson language will become 
d(\) or d(\). Thus we can read off the tachyon operators by looking at modes of p(a;, t\ 
The normal ordering is with respect to the Fermi surface, but this is only relevant 
for the zero energy tachyon which we will not consider. 

It is straightforward to write down the tachyon operators from the normal or- 
dered product of two fermions. We have the operators 

t 



Tout,fe = y rfcu^out,..-fc^L,^- (3.27) 

These have the property that for A; > 0, tI^^. creates a positive energy incoming 
tachyon of momentum fc, whereas Tout,fc acting on creates an out state with 
a positive energy outgoing tachyon of momentum k. Clearly, Tj|^^ = Tin _fc and 
similarly for Tout- The commutation relations are indeed those of a free scalar 

[^in,fc, = [Tout,., Ti,,.,] = k6{k - k'). (3.28) 
Up to a rescaling, these are of course the same as the tachyon modes that we wrote 



m { 

We now have deflned brane creation and annihilation operators, and tachyon 
creation and annihilation operators. We can now write down a general amplitude 
involving branes (with prescribed open string tachyon proflles) and closed string 
tachyons, and the matrix model result for such an amplitude can simply be worked 
out used the commutation relations of the \&in,out operators, and the fact that 

^k^lAi) = *Lt,^lAi) = for u>ii 
= ^out,cu|/u) = for u; < /i. 

Notice that it matters whether we flrst create a brane and then a tachyon, or 
whether we first create a tachyon and then a brane. 
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This is a signal of the fact that there is duahty between bosons and fermions 
(bosonization), or equivalently, that there is a duahty between open and closed 
strings. The complete Hilbert space can be spanned either by closed strings states 
plus their solitons (the closed string = boson picture), or by systems of branes and 
anti-branes (the open string = fermion picture). 

To illustrate this picture, we will now compute a few quantities and check the 
results against results in the literature and our cylinder calculation described in 
section (0). 



Tachyon correlators have been computed in e.g. [17| with precisely these con- 



ventions. The world-sheet tachyons are obtained by including a leg-pole factor, e.g. 

in,aj;world— sheet -p/ ^(^j) \'J-<j^j 



In the conventions of we should include an extra factor of y/27r in the right hand 
side of this expression. One also has the option of including a factor on the right 
hand side, this is a simple choice of normalization, and one could also include phases 
in the definition but we will not do that here. 

The first quantity we will consider is the emission probability of a tachyon by a 
rolling tachyon. For this purpose we need to compute 



(3.31) 



This can easily be computed using Wick contractions, and we find 

r■min{^l,^l+k) ^^^-jo; log Aout-i(a;-fc) log _ 

71 ^ ^ —■ (3.32) 

If we compute numerator and denominator in the stationary phase approximation, 
which is the leading contribution in the l//i expansion as discussed above, we find 
for the emission probability 

_g-*feiogAout^ (3.33) 

Up to a prefactor which depends on the choice of convention, this is the answer 
one expects for the emission probability for a tachyon from the decaying brane. 
Inserting T^y^ would naturally describe the absorption of the tachyon by a decaying 
brane. We could also have chosen the order of the out brane and out tachyon operator 
the other way around. Then the results of the calculation would have no tree level 
contribution. 

A similar calculation yields for the tree-level contribution to the emission ampli- 
tude 

, /^|-^out,Aout-^in,fc-^in,Ai„ ■ 



r~\ out,Aout m,Ai„"^ 



_g-*feiogAi„ ^334^ 
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Inserting Tin again would describe absorption of a tachyon by an incoming brane. 
3.2 Two branes 

Finally, we consider the situation with two D-branes. From the boundary state point 
of view, we have one in-boundary state corresponding to one D-brane, and another 
out-boundary state corresponding to the complex conjugated D-brane. Therefore, 
we propose that the relevant matrix model calculation that describes the interaction 
between two D-branes is 

^=(^out,AL.An.A?^L.. Dlu) (3.35) 



The operators D. \2 and D'^ that describe the second brane have been complex 
conjugated, so that the role of in and out has been interchanged, and their definition 
therefore involves phases e^'^^°^^^^ and e"^'°^^°"t. 

The amplitude in question can be rewritten using Wick's theorem as 

J —oo J —oo 

- r du;e*^('°s^i"-'°sA^„) r (^^e*^('°s^-*-'°sALt) (3.35) 
J —00 J —00 

The first term in the expression has the interpretation as interactions between the 
incoming and outgoing branes, the second term gives the interaction between the 
incoming-incoming and outgoing-outgoing branes. More precisely, in terms of the 
boundary state computation, each integral corresponds to the interaction between the 
relevant incoming/outgoing states. The incoming/outgoing states are described by 
choosing the relevant time contour, or picking the part of the full brane computation. 
Thus the above four-point calculation gives the "square" of the cylinder. Also notice 
that the above result contains full information of the interaction beyond cylinder 
diagram. The information is contained in the reflection coefficient 

To compare with the string theory calculation therefore we take a single integral 
piece. Focusing on the piece in the second term which describes the interaction of 
two incoming branes, we have 



•/ — oc 



J —00 

which can be evaluated to be 

g.M(iogA?„-iogA^j Lsi^iogXt^ - log AL) - , \ . (3.38) 

The presence of the (5-function is understood from the exclusion principle, one cannot 
place two fermions on the same location. Assuming Ai is not equal A2, and neglecting 
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overlap with the Fermi-sea we get for the cyhnder diagram between two incoming 
branes 

} (3.39) 



exp -7e + 







^ ifc(logA2^-logAiJ 
k 



This result has the interpretation as the exponential of a single string theory annulus 
diagram, in agreement with the computations in section ^ Note that the full compu- 
tation of ( p.36|) has the other piece multiplied to (|3.39|) , which has the interpretation 
as the annulus diagram between outgoing states. Thus ( ^.361 ) can be thought of as 
the combination of two annulus diagrams as depicted in Figure |1|. The two point 



, out 



< in 



out 




Figure 1: Interactions of two DO-branes are composed of two terms each of which contain 
two cylinders. The shaded area denotes the strongly coupled region. 

function calculated in section |^ corresponds to the contribution of one of the two 
cylinders. We see that the cylinders are originated only from the asymptotic infini- 
ties. This is because these cylinders are the genus suppressed ones as we explain in 
Appendix A. Closed string modes emitted or absorbed at the strongly coupled region 
will be affected by higher genus corrections. 

Notice that (|3.36|) is the exact answer in the matrix model, and that again we 
see no sign of divergences and/or the exchange of discrete states. 

4. Annulus diagram from boundary states 
4.1 Preliminaries 

In this section, we consider the interaction of two unstable DO-branes in terms of 
boundary state formalism, and compare with the results of the previous sections. 
Each brane carries an open string tachyonic degree of freedom. The rolling of the 
tachyon can be described by boundary conformal field theory with a boundary in- 
teraction 

S = ^j dte'''^^ (4.1) 

^Note that the integration is divergent when k approaches zero. We regulate the divergence, 
and Euler's constant 7e appears by the process of regularization. 
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This interaction describes a decaying (half) brane solution, where the tachyon rolls 
down its potential. One can also consider a bounce (full brane) solution with a cosh 
boundary interaction 

S = X [ dt cosh X%t) (4.2) 



This perturbation can be interpreted as a brane building up from a finely tuned 
incoming coherent tachyon state, and decaying again by rolling down the potential. 

In the matrix model, the decaying solution is interpreted as an eigenvalue rolling 
down the potential, while the bouncing solution as an eigenvalue coming in from 
asymptotic infinity, rolling up the potential, then rolling down again. The precise 
interpretation of these boundary states in terms of free fermions was discussed in 
detail in section |. 

The boundary state in the two dimensional string theory decomposes into an 
dependent (rolling tachyon) part, a Liouville part, and a ghost part 

\B) = \B)xo \B)l |fi)ghost (4.3) 

The unstable DO-branes are described by the ZZ boundary states in the Liouville 
theory [O, 113. The rolling part of the boundary state is less well understood, and 



we first turn to a brief discussion of the rolling part of the boundary state. 
4.2 The boundary state for the rolling tachyon 

We are interested in determining the boundary state that describes the boundary 



interaction ( ^4.1| ) in the free theory that consists of a timelike scalar. This is some- 
what subtle issue. Two alternative representations of this boundary state have been 
proposed in the literature. In the first one we expand the boundary state in Ishibashi 
states \uj)) built out of normalizable operators in the free timelike scalar theory. This 
representation is naturally obtained by considering a spacelike Liouville theory with 
arbitrary background charge, and by taking a limit where the background charge 



vanishes. This leads to the boundary state p2 



\B)xo = / due-'-'^^^'-^lu;)), (4.4) 
J smh iruj 

where 

A = vrA. (4.5) 

There is an alternative representation of the boundary state, obtained by Wick ro- 
tating the boundary state for an Euclidean boson with boundary interaction e*"^; 
this boundary state is an expansion in Ishibashi states built on non-normalizable 



operators and its form can be obtained from |25, 26 



^)^° = EE (^tr) (4-6) 

j m>0 ^ ^ 
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Here the sum is over nonnegative half- integer j, and the sum over m is such that 



m + j is an integer. For some more discussion of these boundary states see e.g. 



and references therein. This representation of the boundary state was used in the 



original work of Sen [15 



Both representations ( ^4.4| ) and ( ^.6| ) are obtained by indirect procedures, either 
taking a limit of Liouville theory or by a suitable analytic continuation. One can 
also try to directly compute the boundary state by considering disk diagrams with 
an operator inserted in the middle of the dis, see e.g. ^ These calculations 
are in principle well-defined, except for the contribution of the zero-mode, which is 
where all the subtleties arise. If we represent the operator inserted in the middle of 
the Dis through the state-operator correspondence as 

for some polynomial V, then the disk calculation will lead to an answer of the form 
(where t is the zero mode of X°) 

/oo 
dte*"*5]]A"e"*c„(P) (4.8) 

71=0 

where Cn(V) does not depend on t or A. It represents the disk calculation with n 
insertions of the boundary operator and a zero momentum operator V{dX^, dX^, . . .) 
inserted in the middle of the disk. We can now see how we have to treat ( [4.8|) in 
order to recover either (|4.4|) or (|4.6|) . 

To recover ( [1.6|) , we should apply momentum conservation to ( [1.8|) , and interpret 
the t integral as 

e^'^*e"*dt ~ 5(cj - zn) (4.9) 



which one would naturally get by Wick rotating the zero mode and viewing the 
integral as a Fourier integral. With this interpretation, there is only a contribution 
for imaginary momenta u = in and this leads to the boundary state ( [4.61 ). The 
coefficients Cn{V) have been computed for various P in f^, leading for example to 
the results 

c„(l) = (-1)", c„(9X°, BX^) ~ (-1)" - 26n,o. (4.10) 

The Snfl that appears is directly related to the nontrivial combination of Ishibashi 
states present in ([4. 6]). 

To recover ([4.4| ), we first consider V = 1. Then (|4.8|) becomes 



with 



j dte''^'p{X,t) (4.11) 



p{X,t) = —^ (4.12) 
1 + Ae* 
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In order to make sense of ( 4.11 ), we have to choose an integration prescription. In 



21[] the integral is defined by taking the contour over the real axis with a small 



rotation t ^ t{l — ie), to include the poles of p{t). With this definition we find 

dte*"V(A, t) = e-^"'°g^ (4.13) 

smh Tcuj 



which is precisely what we have in ( [4.4| ). However, we only considered the case V = 1, 
and we should extend this calculation to include all V. One can see from simple 
examples that this is not at all straightforward. For example, taking V = d'^X^d^X^ 



[|] leads in (g]|) to 

j dte'^^'ipiX, t)-2 + Ae*) (4.14) 

and the contour integral given above does not yet give a good definition of the integral 
J (ite*'^*e*. One way to define these integrals is to demand that the resulting boundary 
state that we obtain is conformally invariant, i.e. obeys 

(L„-Z_„)|5)xo = 0. (4.15) 

Since the highest weight representations of the Virasoro algebra with c = 1 and 
conformal weight —u"^ < do not contain any null vectors^ ^, it is sufficient to find 
the disk diagrams with V = 1 only; all other disk diagrams then follow from ( |4.15| ). 
Thus, conformal invariance of the boundary state plus the prescription ( [4.13| ) provide 
a unique answer for all ill-defined integrals such as ( [4.14| ), and the resulting answer 



is then given by (|4.4|) . All this is true except at zero momentum. For uj = 0, 
the corresponding highest weight representation has null vectors and (|4.15|) is not 
sufficient to determine the boundary state. This will not play an important role in 
what follows, but we note that the contribution to ( |4.6| ) with zero momentum is the 
sum of the terms with m = 0, so it simply is \ j, 0, 0)), which behaves as if there 
were no null vectors anyway. 

We have now motivated the two representations ( [4.4| ) and ( [4.6| ) from direct world- 
sheet calculations. It appears that some of the information of ([4. 61) got lost in 



relating it to ( [4. 41) ; in particular, we could have taken the zero mode approximation 



of ( [4. 61) as in [|T5|, and from that we could also have derived (|4.4| ) via a Fourier 



transformation. On the other hand, the zero mode structure of ([4.6|) plus conformal 



invariance determine (4^) up to the contribution of = 0, so perhaps it should be 



possible to somehow re-obtain ([4.6|) from ( [4.4|) as well. 



In this paper we are mainly interested in computing an annulus diagram involving 
two boundary states. Though ([4.4|) and (WM) are related, it is not clear whether and 



^The Kac determinant for c = 1 theories at level N is Aj^ih, c) - npg<A'('*- (p-'?)V4)^^^~^'^ 
and clearly this has no zeroes for h < 0. 

^The absence of null-vectors is related to the statement in |Q that one can choose a gauge at 
nonzero energy where the timelike oscillators are not present. 
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how calculations involving (|4.4| ) are identical to calculations involving (f4.6|) . If we use 
( [4 .61) one might be inclined to write down the following expression for the annulus 
amplitude 

j m>0 ^ ^ 

It is not clear that this is correct, since we assumed here that the Ishibashi states 
|j, m, m)) are orthonormal to each other. The zero- mode contribution to the overlap 
of two such Ishibashi states involves integrals of the form J dte*^"""''-'* which is not 
well-defined. Taking it equal to 5n,n' leads to (|4.16|) . However, it is not obvious that 
this the correct thing to do. Another peculiar feature of ( [4.16| ) is that it contains 
binomial coefficients that grow very rapidly. All these problems are related to the 
fact that ( [4.6| ) is an expansion in non-normalizable states, and a general framework 
to formulate string amplitudes in terms of such states does not yet exist. 

The boundary state ( [4. 41 ) on the other hand is a much more conventional expres- 
sion in terms of normalizable Ishibashi states. It is this type of expansion that we 
would normally use when computing string theory amplitudes, and we will use (|4.4| ) 
in the rest of this section. It turns out that ([4.4]) yields answers that can be directly 
compared to matrix model calculations. It remains to be seen whether and how the 
same results can be recovered from ( |4.6|) . 

Except for the half-brane (^4.1|) we will also be interested in the full-brane ( [4.2| ). 
For the full brane a discussion similar to the one above applies, the main difference 
being that p(A, t) in ( [4.12| ) gets replaced by 



p(A, t) = — K— + \ 1, A = sin Att. (4.17) 

1 + Ae* 1-Ae-* 



For the full brane we need to choose in ( 4.13|) an integration contour. One can 



choose a real contour, which corresponds to a brane building up from a finely tuned 
tachyon state, and then decaying back. Another contour choice is the Hartle-Hawking 
contour, consisting of a piece along the imaginary time axis from t G (0, ioo), and a 
real part t G (0, oo) (or t G (— oo,0)). These correspond to a preparing a state at 
t = and evolving it to (from) the asymptotic infinity pT[] . 

4.3 Cylinder computation 

We now concentrate on computing the cylinder diagram of two decaying DO-branes. 
In terms of boundary states it can be written as 

Azz = \l — (fel® (^A'lg'^'I^A)® |5zz)Aghost (4.18) 

In this expression \B\) is the time dependent part of the boundary state, and \Bzz) 
is the Liouville part. The ghost contribution is separated into Aghost- We consider 
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the overlap of two different rolling boundary states, parametrized by A and A'. Also, 
q is the modulus of the annulus. Since the interacting states are time-dependent, we 
compute the full time-integrated diagram. 

The Liouville part of the annulus diagram is already given in [14, 23], 

Zi,i = J dpxp^i,i{p)^i,i{-p) (4.19) 

where the index {1, 1} refers to the fact that we take the m = 1, n = 1 boundary 
state, \E'i^i is the wave function, and the character is 



Xp 



qP 



r]{T) 

Vir) =q^^l[il-qn (4.20) 



oo 



n=l 



More explicitly 
The ghost part of the annulus is 



A^M. = dMf (4.22) 



We still need to compute the time- dependent part. Inserting Ishibashi states built 
over e^^^^lO) 

Ax,y = (5v|g2^°|5A) = jdtj dt' j du j d^'(5v|^)((^|g2^'V'))(^'|5A) 

(4.23) 



.,2 



du^piX, t)p(X', t')e*"(*-*') (4.24) 

v{q) 

As we discussed in the previous subsection, we used for our boundary state 
one that is expanded purely in terms of normalizable Ishibashi states, and which 
is therefore determined by the zero-mode structure of the corresponding Euclidean 
boundary state, together with conformal invariance. 

Assembling the full annulus amplitude and performing the integration over the 
modulus q we obtain 

A.. ^-Ijitj M'j ± j ,.p(A.<)p(A',0^5;«e""-> (4.25) 

Note that the ghost contribution canceled the //-functions. We now have a choice 
to proceed with the u or with the t integration first. Let us proceed with the u 
integration, 

^zz = - / * / *' / f P(A. t)p{y, t) (2l!^) ' (4.26) 
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Finally for the time integrations the contour has to be specified. The integral of 
p{t) has poles at discrete values on the imaginary axis. For the half-brane we follow 
the prescription described in the previous section, with the result 



dtp(t)e'P' = e-'P^°^^ (4.27) 

smh np 



So finally we arrive at 

Azz= I 

P 

This result is in agreement with the matrix model computations p.l3| and |3.39 . 



A^^ = [ ^e-'PiiogX-iogX') ^4 28) 



If we consider fuU-brane or bounce solution with real contour, the annulus am- 
plitude is modified to 



— J ^ ^g-«plogA _ g+jplogA^ ^g+iplogA"' _ g~iplogA'^ 

Had we done the t integration first, we would have obtained instead 



(4.29) 



Here we have taken the half-brane case for simplicity. Since the p integral diverges, 
we now have to proceed with the uj integration. The expression has double poles 
at discrete imaginary values of uo, which indicates the presence of discrete modes. 
Performing the uj integration we get the previous result ( 4.28| ) plus an infinite sum 



dp 1 f sinhvrp 



containing the contribution from the discrete modes. This infinite sum with contri- 
butions from discrete modes does not seem to be refiected in the matrix model result. 
The expression ( [4.28 ) contains only an exchange of intermediate tachyons, and this 



result is in agreement with the matrix model computations. 



5. Discussion 

We have seen that the matrix model computations performed in the free fermionic 
picture agree with the string theory cylinder diagram in which only the tachyon 
exchange is taken into account. The matrix model correlators in the tree level {qs 
0) approximation do not indicate the presence of discrete modes. 

On the other hand, one expects a sign of the presence of discrete states, at least 
in Euclidean signature, where they correspond to normalizable string states. In that 
case they are part of the full spectrum, and carry important physical information 
such as gravitational interactions. The discrete modes certainly can be seen to be 
implicitly present in the matrix model through the ground ring generators and the 
Woo algebra [E^ . 
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In the cylinder computation between the two DO-branes we have encountered 
several difficulties. A better understanding of these difficulties might provide inter- 
esting insights into the nature of the matrix model, in particular regarding the role 
of the discrete modes and its non-perturbative degrees of freedom. In particular, we 
would like to know whether the cylinder calculation can be reformulated purely in 
terms of ( [4. 61 ) or not. Another confusion arises relating a Euclidean worldsheet to a 
Minkowski spacetime computation. 

Another inherent difficulty with the cylinder computation is the divergence of 
the p-integration. We have chosen to take integrals in a particular order to avoid 
the divergences. This corresponds to regulating the divergent integral in a specific 
way. We have chosen the regularization to be in agreement with the matrix model 
computation, but again, this regularization procedure cuts off an infinite sum of 
discrete modes. Clearly, a better understanding of the string theory computation is 
desirable. 

There are several ways in which the results in this paper can be extended. First of 
all, we can consider the interactions of p > 2 branes. From the matrix model point of 
view, this can also be understood as taking (N+p) x (N+p) matrices, integrating out 
the NxN degrees of freedom and taking a large limit. According to [19|, the result 
is given by a Kontsevich matrix integral whose result is a certain r-function. This 
r-function has an alternative representation as a determinant of fermion bilinears, 
which is precisely what one obtains from a straightforward generalization of ( 3.35| ). 
This provides further evidence for the structure advocated in section |^. It would 
be interesting to understand more directly how such a structure can emerge from 
world- sheet calculations. 

It is also relatively straightforward to generalize our results to rolling tachyons 
in the c = 1 model of [|l^ . This simply amounts to filling both sides of the potential, 
and incorporating both asymptotic regions of the matrix model. In particular, we 
can study rolling tachyons corresponding to perturbations of the form AsinhX°, 
which start at early times at x — >• — cxo and at late times end up at x — > +oo. It 
is easy to write down corresponding expressions in the matrix model, but they are 
more difficult to understand on the world-sheet, as these are configurations that 
penetrate deeply in the non-perturbative regime; they start out as coherent states of 
the symmetric combination of the closed string tachyon and RR scalar, and end up 
as coherent states of the antisymmetric combination instead. 

Our calculations also clearly support Sen's point of view [2^ that open and closed 
strings are dual pictures of one and the same thing. The matrix model as expressed in 
terms of the original variables (eigenvalues) is the open string picture. The fermionic 
description, which is a reformulation of the original matrix model, represents another 
formulation. The dual closed string picture arises from bosonization of the fermions. 
In the closed string picture there is a clear distinction between closed string states 
and D-branes: the second are coherent states of the ffist (in other words, they are 
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solitons in the closed string). Similarly, in the open string picture closed strings 
arise as excitations of brane/antibrane systems with total brane number zero. Thus, 
both in the closed as well as in the open picture, there is a conserved fermion=brane 
number, which is conserved despite the fact that the branes are unstable. 
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Appendix A. Stationary phase limit of the two point function 

In this appendix, we consider the full two point correlator of the free fermions with- 
out taking any asymptotic bosonization. At the genus zero limit, we can take the 
stationary phase approximation. This approximation extracts the piece of the corre- 
lator corresponding to the tree level {gg 0) string computation discussed in section 



The full fermionic two point correlator is |18| 

(/i|r(V^(ti, Xi)V^(t2, = J dio {e{uj - fi)9{At) - 9{fi - iu)9{-At)) e""^* 

X {i;'^{uj,Xi)ip^{uj,X2) +i'^iuj,Xi)i)~{uj,X2)) (5.1) 

The subscript E refers to taking the Euclidean two point function. Using the proper- 
ties of wave functions of the inverse harmonic oscillator, or by analytically continuing 
the harmonic oscillator propagator we arrive at the equivalent expression ||18 



^„ /•sgn{p)-oo sp+ifis 

(/i|T(^(ti,xO^(t2,X2)|/x)^ = te-^""' / ^e-'''^' / ds- 



271 Jq (-47risinhs)i/2) 
/ i ((x? + Xo) cosh s - 2x1X2) \ ,^ 
" [ I^s ) ^'-'^ 

To consider tree level interaction, we suppress the higher genus corrections. We 
rescale x /5^/^x, and /i Pfi, where P is sent to infinity and Pjj is kept fixed 
but large. Taking this limit suppresses both loop corrections (these are effectively 
h = 1/(3 corrections), and tree level (l//5/i) higher order interactions. The first 
corrections would correspond to the higher genus corrections in the string diagrams. 
The tree level corrections would come from taking into account the interactions of 
the bosonized scalar field. Since in this paper we are mostly interested in the cylinder 
computation on the string side, we keep Pfi ^ 1. 

Rescaling x the integral acquires a large phase, and can be evaluated with the 
stationary phase method. The phase will be stationary at 

e^- + e-^=^ + ^ (5.3) 

X2 Xi 
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where s denotes the stationary value of s. Depending on the signature of p in the 
integration and on the signature of (xi — X2), this gives 

s = Q{p) ( e{xi - X2) log (— ) + e(a;2 - xi) log ( — ) 

\ X2 Xi 

+ 0(-p) f 0(a;i - X2) log (^) + 0(X2 - xi) log i^)) (5.4) 

Here G is the step function. At the stationary point part of the Gaussian factor 
coming from the integration of the quadratic part cancels the 1 / V sinh s factor in the 
integral. The rest of the expression assembles to the WKB wave function, and an 
additional exponential factor, giving 



s{t,, t2, xi, X2)^ = ixhr^^^'^'""' r ^e-^^^* {e(p) [e(xi - x^) x 

y p\xiX2\ Jo 27r 



+ Q{-p) 



6(xi-X2)e 41.^2 ^iJg^^f"^'^ 



+ e(x2-Xi)e-x(-?-i)e(*''^-^)^°sSj| (5.5) 

Going back to Minkowski signature requires At iAt, p ip. The four pieces of 
the integral correspond to picking the left or right moving pieces in the correlator: 
R^L,R- Dropping the WKB factor and choosing the left-left piece we obtain for 
the two point correlator 

S{h,t2,x„X2)= r dpeM'^^^.-^') (5.6) 



Taking the orbits Xi^2 = Ai^2Cosht (where t is sent to infinity) we obtain the 
result (|2.13| ) computed from the asymptotic bosonization. 

Notice that to get this result, we have only suppressed higher genus and higher 
order tree level correction, we did not consciously remove any possible discrete modes. 



Appendix B. Two point function of eigenvalue density 

In this appendix we consider two point function of the eigenvalue density. Roughly 
speaking, since for chiral fermions ipi^x^t) ~ exp (J ip'^ipi^x^t)) , by computing the two 
point function of the eigenvalue density we can try to see once more whether discrete 
modes are exchanged between two D-branes, and we can also present an alternative 
derivation of the cylinder diagram 

The two point function of eigenvalue density is defined in terms of 

G{ti,Xi]t2,X2) = {fJ^\lp''i^{ti,Xi)lp^ 1^(1:2, X2)\fi) 

G(cui,Xi;c^2,X2) = [ dtidt2e''^'''e''^'''G{t,,Xi;t2,X2). (5.7) 
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We can compute this using the two point functions of fermions computed in ap- 
pendix A, 



'0 -'-oo S + t 

2isH\^ \ I ^ \JlitH\ 



(xi|e^*^^|x2)(x2|e^'*^|xi), (5. 
where {xm\^^^^\xri) defined as 



^ ™^ ^ (— 47ri sinh s)-*^/^) ^\ 4sinhs / ^ ^ 

Here we have taken the Euchdean time. As in the Appendix A, We will think 
of the genus zero limit, scaling /i — \f^x and taking the limit /? — oo 
with (3^ = fixed ^ 1. Then we can apply the stationary phase approximation 
to the s and t integration. But as before, there are two sets of saddle points : 
s = i = ±log(xi/x2). To have a non- vanishing denominator in ( |5.(j| ), we choose the 
same sign for the saddle of s and t. Finally, we get 

ZX1X2 log(a;2/a;i) 

We take the inverse Fourier transformation to (t, x) coordinate space 

G{ti,xi-t2.X2) = j duJiduo2e-'^'''e-'''^'''G{uJi,xi-U2,X2) (5.11) 
Stripping off WKB wavefunctions, and Wick rotating to Minkowski space, we 

get 

duj ^ g-i<^(<i-logxi-t2+logX2)_ ^5^^2) 

logXi — logX2 

To get the cylinder diagram between two D-branes, we should integrate over t 
in ( |5.12| ) up to a large nonzero value of t, and then insert the classical orbit at the 
boundary value of t, which yields the simple answer 

Jo ^ 

Since ip{xi,ti)-ip{x2,t2) ~ exp (J ijj^tjj{xx,ti)) exp (J tjj''ip{x2,t2)),hj summing the all 
the ladders, we recover ( |2.13| ). Note that we again see no divergences due to discrete 
mode exchange (or any other sign of discrete modes for that matter). 

We can evaluate the two point function of the eigenvalue density in another way. 
If we take the Fourier transformation from x to z, and rotate z to il, we get the two 
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point function of macroscopic loop operators of loop length /. The result is already 
well-known in the literature [l^ . At genus zero the result is 

We can now analytically continue back to 2; = —il and do an inverse Fourier trans- 
formation to get the following expression in terms of x, 

f rfcuicuie'^-i^*^-*^) X 



cosh 2071 sinh — cosh 2co'i sinh — . (5.15) 



To get the Minkowski result we perform a Wick rotation t —it and u — > iu. Again 
taking off the WKB factors, we have 

~ j dcuicuie-*^i(*^-*^-'°*^^^+^°s^^). (5.16) 

Also here we should do the integral over t to reproduce the correct answer for the 
cylinder diagram, 

A = dtidt2j rfcJiCJie-*^^(*i-*2-^°s^i+^°s^2) 

UJ 

In the last line, we have again inserted the classical trajectory at asymptotic infinity. 
We find that the result agrees completely with the previous results. 
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